Pulling an adsorbed polymer chain off a solid surface 
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The thermally assisted detachment of a self-avoiding polymer chain from an adhesive surface by 
an external force applied to one of the chain ends is investigated. We perform our study in the "fixed 
height" statistical ensemble where one measures the fluctuating force, exerted by the chain on the 
last monomer when a chain end is kept fixed at height h over the solid plane at different adsorption 
- - - strength e. The phase diagram in the h — e plane is derived both analytically and by Monte Carlo 

' simulations. We demonstrate that in the vicinity of the polymer desorption transition a number of 

, properties like fluctuations and probability distribution of various quantities behave differently, if h 

' rather than / is used as an independent control parameter. 
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^ '. I. INTRODUCTION 

The properties of single polymer chains at surfaces have received considerable attention in recent years. Much of 
I this has been spurred by new experimental techniques such as atomic force microscopy (AFM) and optical/magnetic 
^ ' tweezers [p which allow one to manipulate single polymer chains. Study of single polymer molecules at surfaces, 
• ' such as mica or self-assembled monolayers, by Atomic Force Microscopy (AFM) method provides a great scope for 
. experimentation 0, i, i, i, i, 0, i- Applications range from sequential unfolding of collapsed biopolymers over 
Ch ' stretching of coiled synthetic polymers to breaking individual covalent bonds[9l.[l0l.[l]|. 

I , In these experiments it is customary to anchor a polymer molecule with one end to the substrate whereas the other 
' end is fixed on the AFM cantilever. The polymer molecule can be adsorbed on the substrate while the cantilever 
' recedes from the substrate. In so doing one can prescribe the acting force in AFM experiment whereas the distance 
between the tip and the surface is measured. Conversely, it is also possible to fix the distance and measure the 
corresponding force, a method which is actually more typical in AFM-experiments. From the standpoint of statistical 
mechanics these two cases could be qualified as /-ensemble (force is fixed while fluctuating the height of the chain 
end is measured) and as ^.-ensemble [h is fixed while one measures the fluctuating /). Recently these two ways of 
descriptions as well as their interrelation were discussed for the case of a phantom polymer chain by Skvortsov et al. 
p^ . In recent papers J3, 14] we studied extensively the desorption of a single tethered self-avoiding polymer from a 
flat solid substrate with an external force applied to a free chain's end in the /-ensemble. 

In the present paper we consider the detachment process of a single self-avoiding polymer chain, keeping the 
distance h between the free chain's end and the substrate as the control parameter. We derive analytical results 
for the main observables which characterize the detachment process. The mean value as well as the probability 
' distribution function (PDF) of the order parameter are presented in close analytical expressions using the Grand 
, Canonical Ensemble (GCE) method p^. The basic force-height relationship which describes the process of polymer 
ILjl ' detachments by pulling, that is, the relevant equation of state for this system is also derived both for extendible as 
well as for rigid bonds and shown to comply well with our Monte Carlo simulation results. We demonstrate also that 
a number of properties behave differently in the vicinity of the phase transition, regarding which of the two equivalent 
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H ' ensembles is used as a basis for the study of systems's behavior. 
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II. SINGLE CHAIN ADSORPTION: USING DISTANCE AS A CONTROL PARAMETER 



A. Deformation of a tethered chain 



Before considering the adsorption-desorption behavior of a polymer in terms of the chain end distance h, we first 
examine how a chain tethered to a solid surface responds to stretching. This problem amounts to finding the chain 
free end probability distribution function (PDF) PN{h) where N is the chain length, i.e., the number of beads. The 
partition function of such a chain at fixed distance h of the chain end from the anchoring plane is given by 

Et,niN,h)^t^loPNih) (1) 
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where (3 = 1 — 71 and the exponent 71 = 0.680^15). Here /Z3 is a model dependent connective constant (see e.g. 
ref.[l^). In Eq. ^ Iq denotes a short-range characteristic length which depends on the chain model. Below we 
discuss the chain deformation within two models: the bead-spring (BS) model for elastic bonds and the freely jointed 
bond vectors (FJBV) model where the bonds between adjacent beads are considered rigid. 



1. Bead-spring model 

The form of Pn (h) has been discussed earlier [l^, |T3| and later used in studies of the monomer density in polymer 
brushes p^ . Here we outline this in a way which is appropriate for our purposes. The average end-to-end chain 
distance scales as Rn = IqN" , where is the mean distance between two successive beads on a chain and v « 0.588 
is the Flory exponent [l^. The short distance behavior, h Rn, is given by 
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where the exponent ( « 0.8. For the long distance behavior, h/ R^ ^ 1, we assume, following Ref.fT8|, that the PDF 
of the end-to-end vector r is given by des Cloizeaux's expression j2(| for a chain in the bulk : (r) = {1 / Rn) F {r / Rn) 



where the scaling function F{x) oc x* exp[—Dx ], and the exponents t — {13 — d/2 + vd)l(l 



5^ 1/(1-1/). Here 



and below d denotes the space dimensionality. One should emphasize that the presence of a surface is manifested only 
by the replacement of the universal exponent 7 with another universal exponent 71 (as compared to the pure bulk 
case!). By integration of Pn{y) over the x and y coordinates while h is measured along the z— coordinate, one obtains 
PN{h) oc {h/ Rn)'^^*~^ exp[— I?(/i/i?jv)*]- As the long distance behavior is dominated mainly by the exponential 
function while the short distance regime is described by Eq. ([2]), we can approximate the overall behavior as 
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A comparison of the distribution, Eq. ([3]), with our simulation data is shown in Fig. [T] The constants A and D in 
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FIG. 1: Probability distribution PNih) of chain end positions h above the grafting plane for a polymer with = 128 monomers 
at zero strength of the adsorption potential e = 0.0. In the inset the MC data for PN{h) (solid black line) is compared to the 
theoretic result, Eq. Q. Dashed line denotes the expected slope of ^ « 0.78 of the probability distribution for small heights. 



Eq. ^ can be found from the conditions: J PN(h)dh = 1 and J h^PN{h)dh — i?^. This leads to: 
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where 6 w 2.43 and C ~ 0.8. One gets thus the estimates A « 2.029 and D « 0.670. 

The free energy of the tethered chain with a fixed distance h takes on the form Ftaii(Af, ^) — — fesTlnStaiK-ZV, /i) 
where denotes the Boltzmann constant. By making use of Eqs. ([T]) and ([3]) the expression for the force Jn, acting 
on the end-monomer when kept at distance h is given by 
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One should note that at h/Rfq ^ 1 we have h oc Rn{Rn Jn / ksTY^ '^^^^^ which, after taking into account that 
= 1 — 1^, leads to the well known Pincus deformation law: h cx loN{lafN/kBTy^'^~^ |2l|. Within the framework of 
this approximation the (dimensionlcss) elastic energy reads Uci/ksT — —N{lofN /kBTY'^ . In result the corresponding 
free energy of the chain tail is given by 
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Eq. ([6]) indicates that there exists a height /iq = [Q/ 50)^/^ R^ over the surface where the force changes sign and 
becomes negative (that is, an entropic repulsion dominates). According to Eq. ^ the force diverges as /jv oc —ksT/h 
upon further decrease of the distance h. 



2. Freely jointed chain 

It is well known pj| that the Pincus law, Eq. ^ , describes the deformation of a linear chain at intermediate force 
strength, l/N"^ ^ lofN/ksT < 1. Direct Monte Carlo simulation results indicate that, depending on the model, 
deviations from Pincus law emerge at h/R^ > 3 (bead-spring off-lattice model) or H/Rn > 6 (Bond Fluctuation 
Model [2^). In such "overstretched" regime (when the chain is stretched close to its contour length) one should 
take into account that the chain bonds cannot exparrd indefinitely. This case could be treated, therefore, within the 
simple freely jointed bond vectors (FJBV) model [22, where the bond length is fixed. In this model the force - 
deformation relationship is given by 

ksT 1 / h 



where denotes the inverse Langevin function C{x) — coth(x) — 1/x and Iq is the fixed bond vector length. 
We discuss the main results pertaining to the FJBV model in Appendix A. The elastic deformation energy reads 
Uci/ksT = —{lofN/ksT) J2iLi < cos9i >= -~N{lofN/kBT)C{lofN/kBT), where 9i is the average polar angle of the 
«-th bond vector (see Appendix A). Thus the corresponding free energy of the chain tail for the FJBV model reads 

f^ = -iVe(^)-^lnM3 (9) 

where we have used the notation G{x) — xC{x) — 2;coth(x) — 1. Now we are in a position to discuss the pulling of 
the adsorbed chain controlled by the chain height h. 



B. Pulling controlled by the chain end position 



Consider now an adsorbed chain when the adsorption energy per monomer is sufficiently large, e > Sc, where 
Ec denotes a corresponding critical energy of adsorption. Below we will also use the notation e = e/kBT for the 
dimensionless adsorption energy. The problem of force-induced polymer desorption could be posed as follows: how is 
the process of polymer detachment governed by the chain end position hi Figure (^ti- gives a schematic representation 
of such a system, and the situation in a computer experiment, as shown in the snapshot Fig. [^b, is very similar. 

As is evident from Fig. [2^, the system is built up from a tail of length M and an adsorbed portion of length N — M . 
The adsorbed part can be treated within the GCE approach [3] . In our earlier treatment [3| it was shown that the 
free energy of the adsorbed portion is -Fads = ksT^N — Af ) lnz*(e), where the fugacity per adsorbed monomer z*(e) 
depends on e and can be found from the basic equation 



(10) 
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FIG. 2: (a) Schematic graph of an adsorbed polymer chain, partially detached from the plane by an external force which keeps 
the last monomer at height h. The total chain is built up from a tail of length M and an adsorbed part of length A'^ — M. 
The force /ai acting on the chain end is conjugated to h, i.e., /m ~ dFtaii/dh. (b) A snapshot from the MC simulation: 
iV = 128, h = 25.0, e = 4.0 and {/) = 6.126. 



The so called polylog function in Eq. PH)) is defined as z) = X^J^i -^"/'^" ^"^^ connective constants /xs, 
^2 in three and two dimensional space have values which are model dependent [T5| . The exponents a = \ + (f) 
and A = 1 — 7^=2 where (p « 0.5 is the crossover exponent which governs the polymer adsorption at criticality, 
and in particular, the fraction of adsorbed monomers at the critical adsorption point (CAP) e = Ec. The constant 
ld=2 = 1.343 Finally w = exp(e) is the additional statistical weight gained by each adsorbed segment. 

In equilibrium, the force conjugated to h, that is, /m — dFtnn/dh, should be equal to the chain resistance force to 
pulling fp = (JtBT /lQ)J-{e) (where J-{t) is a scaling function depending only on e), i.e., 
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¥^^"H^)=^P , for FJBV-model 

The resisting force /p holds the last adsorbed monomer on the adhesive plane (see again Figl2^ whereby this monomer 
is shown to experience a force /a/). One should emphasize that the force /p stays constant in the course of the pulling 
process as long as one monomer, at least, is adsorbed on the surface. Thus /p corresponds to a plateau on the 
deformation curve (force / vs. chain end position h). The adsorbed monomer (see Fig. [5]) has a chemical potential, 
Mads = Inz*, which in equilibrium should be equal to the chemical potential of a desorbed monomer in the tail, 
jjLdes = d{Ftaii/ ksT) /dN . The expression for Ftaii depends on the model and is given either by Eq.® for the BS- 
model or by Eq.® in the case of FJBV-model. Taking this into account the condition /iads = /^dcs leads to the 
following "plateau law" relationship 
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{\ln[fi3Z*{e)]\) , for FJVB-model 



where Q~^{x) stands for the inversion of the function G{x) — a;coth(a;) — 1. One should note that Eq.(fT2| coincides 
with Eq.(3.16) in Ref. [3| which determines the detachment line in the pulling process controlled by the applied 
force. Close to the critical point Cc, the plateau force fp goes to zero. Indeed, since in the vicinity of the critical 
point ln[^3Z*(e)] oc -(e - ecY^"^ (see ref.[14]), and Q~^ix) « (3x)^/^, one may conclude that /p oc (e - ecT^'^ for the 
BS-model and /p oc (e - ej^/^"^ for the FJVB-model. 

One can solve Ea. ([TT|) with respect to M (taking into account that h ^ Rm), and arrive at an expression for the 
tail length 

{if^y^^ ' , for BS-model 
M{h,e) = { (13) 
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where the force at the plateau, /p, is described by Eq. (I12p . If for the degree of adsorption one uses as an order 
parameter the fraction of chain contacts with the plane, n = N^/N , where is the number of monomers on the 
surface, one can write ^J] 
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(14) 



where Fads and Ftaii are free energies of the adsorbed and desorbed portions of the chain respectively. The free energy 
-Fads = kBT[N — M{h,e)]\nz*{e) whereas Ftaii = kBTfidcsM{h,e) (recall that /ides is the chemical potential of a 
desorbed monomer). After substitution of these expressions in eq. (|14p and taking into account that in equilibrium 
Mads = A*dos (the sequence of operations is important: taking the derivative with respect to e is to be followed by the 
condition /iads ~ Mdos) so one gets 
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i.e. the order parameter n is defined by the product of monomer fraction in the adsorbed portion, 1 — M/N, and the 
fraction of surface contacts in this portion, —d\n.z*/de. The expressions for the order parameter can be recast in the 
form 
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Here ci and C2 are some constants of the order of unity. 

As one can see from Eq. (fT6)) . the order parameter decreases linearly and steadily with h/N. This behavior is 
qualitatively different from the abrupt jump of n when the pulling force / is changed as a control parameter. In 
Section |V] we will show that this predictions is in a good agreement with our MC - findings. The transition point on 
the n vs. h curve corresponds to total detachment, n = 0. The corresponding distance h will be termed "detachment 
height" hij. The dependence of ho on the adsorption energy e can be obtained from Eq. (|16p where n is set to zero, 
i.e. 




BS-model 
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(17) 



where again /p as a function of e is given by Eq. (fT2)) . The line given by Eq. ([17]), is named "detachment line". 
It corresponds to an adsorption - desorption polymer transition which appears as of second order since this order 
parameter n goes to zero continuously as h increases. One should emphasize, however, that this "detachment" 
transition has the same nature as the force-induced desorption transition in the /-ensemble where the pulling 
force /, rather than the distance h, is fixed and used as a control parameter. This phase transformation is known to 
be of first order. 

It is easy to understand (cf. with Eg. p^ ) that the condition M{h, e) ~ N corresponds to the detachment line as 
well as to a terminal point of the force plateau. It can be seen in the MC - simulation results in SectiorfVl Fig. [51 



III. PROBABILITY DISTRIBUTION P{K) OF THE NUMBER OF ADSORBED MONOMERS 

The grand canonical ensemble (GCE) method, which has been used in our recent paper fl3|, is a good starting 
point to calculate the probability distribution function P{K) of the adsorbed monomers number K . According to 
this approach, the GCE-partition function of an adsorbed chain has the form 

\ V" " N K Voiwz) Q{z) 
^iz,w)=^ ^^^,KZ w = _ (18) 

N=l K=0 ^ J \ J 



where z and w = exp(e) are the fugacities conjugated to chain length N and to the number of adsorbed monomers 
K, respectively. In Eg. p^)) U{z), V{wz) and Q{z) denote the GCE partition functions for loops, trains and tails. 
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respectively. The building block adjacent to the tethered chain end corresponds to Vq{wz) = 1 + V{wz). It has been 
shown 1J| that the functions U{z), V{'wz) and Q{z) can be expressed in terms of polylog functions, defined in the 
paragraph after Eq. (fTO|. as U{z) — $(a,^3z), V{wz) — <i>(l — 7^=2,^2^2;) and Q{z) = 1 + $(1 — 7i,/i3z), where 
/i3 and fj,2 are 3d- and 2d - connective constants respectively. By making use of the inverse Laplace transformation 
of S(z, with respect to z (see, e.g. [1^) the (canonical with respect to the chain length N) partition function is 
obtained as 



Sjv(i(;) = '^N,K — exp[— In z*(w)] 



(19) 



K=0 



where z*(w) is a simple pole of S(z, w) in complex z-plane given by equation V{wz*)U{z*) — 1, i.e. by Eg. ljlOp . 

The (non-normalized) probability for the chain to have K adsorbed monomers is P{K) oc 'B.n,k = 
exp[—T{K)/kBT], where J'{K) is the free energy at given K. It is convenient to redefine the fugacity it; as w ^ ^lu, 
(as well as Wc ^Wc) where ^ is an arbitrary complex variable. Then the probability P{K) can be found as the 
coefficient of in the function 'E.]\[{£^w) expansion in powers of ^. Therefore 



P{K) = exp[~T{K)/kBT] 
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where the contour of integration is a closed path in the complex ^ plane around ^ = 0. (see e.g. [25|). To estimate 
the integral in Eq. (|20p we use the steepest descent method I25|. 

For large N the main contribution to the integral in Eq. ((20|l is given by the saddle point C = Co of the integrand 
which is defined by the extremum of the function g(^) = — lnz*(^ w) — [{K + l)/N] Inf, i.e., by the condition 



K + 1 



N 



d\nz* 



(21) 



The integral is dominated by the term exp[— In 2:*(^ow) — {K + l)ln^o]- Another contribution comes from the 
integration along the steepest descent line. As a result one obtains 



P{K) oc 
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The validity of the steepest descent method is ensured by the condition of the second derivative Ng"{^Q) being large, 
which yields 





' dliyz*' 
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(23) 



A more explicit calculation whithin this method can be performed in the vicinity of the critical point e = £c 



A. PDF of the number of adsorbed monomers close to the critical point of adsorption 

In this case the explicit form of Inz* is known [14] and after the redefinition of the fugacity, w — > ^w, it reads 

\nz*{^w) = ~ai{w ~ Wc)^/^^^/^ - In^a (24) 

where ai is a constant of the order of unity. The critical adsorption fugacity Wc = exp(ec) is defined by the equation 

C(a) $(1 - 7d=2, Ai2Wc/M3) = 1 (25) 

with (^(x) denoting the Ricmann zcta-function. 

By using Eq. ((24|) in Eg. pTjl one arrives at the expression for the saddle point 
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where 02 = {4>/ai)'^. Using Eq. pi]) and Eq. (l26|l in Eq. ([22|l then yields the expression for PDF 

-4>K-l/2 , w 



PiK)^[iw-w.)e'^^f^-j ^ (27) 

For reasonably large K and after normalization one arrives at the final expression for the PDF 

where we have introduced the usual adsorption scaling variable rj — bi{w ~ Wc)N'^ (61 is a constant of the order of 
unity; see e.g. |26i]) as well as the normalization constant C: 



N 



Civ) -Y^J^K (29) 

K=l 

One can readily see that the width of the distribution increases with w or with e. To this end one may directly 
calculate the fluctuation variance as follows: 

{K - KY ^ -iV ^'i"^' y (X Nw{w/(j) - Wc){w - We)'/^"' (30) 
a(ln w)"^ 

where the expression for h\z*{w) given by Eq. ([M|) (where also <^ = 1 ) has been used. Taking into account that 
4> ~ 0.5, it becomes clear that the variance really grows with w. 

The validity of the steepest descent method is ensured by the condition Eq. (P5|) . Using Eq. ([M|) and Eq. (I26p . one can 
verify that this criterion holds when N{w - WcY{K/NY^'^'>' > 1. We recaU that (p « 0.5, so that {w - Wc)A^^/^ > 1. 
In result the criterion becomes 

^ (w - w,) < 1 (31) 

In the deep adsorption regime this condition might be violated. Nevertheless, the steepest descent method could still 
be used there, provided that the appropriate solution for z*{w) (see Ea. (fTII)) ) is chosen. 



B. The regime of deep adsorption 



In the deep adsorption regime one should use the the solution for z*{w) which was also discussed in ref. [3]. 
Namely, in this case 
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With Eq. ([5^ the mean value K = —Ndlnz*{w)/d\nw can be written as 

1/(1-A) 
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(32) 



(33) 



thus K tends to N with w growing as it should be. The variance of the fluctuations within the GCE then becomes 



{K - KY = -N 



\ii z*(w) 
9(ln w)2 
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1/(1-A) 



(34) 



i.e., the fluctuations decrease when the adsorption energy grows. Comparison of this result with the result given by 
Eq. (|30p leads to the important conclusion that the fluctuations of the number of adsorbed monomers flrst grow with 
e, attain a maximum, and flnally decrease with increasing surface adhesion e. The position of the maximum reflects 
the presence of finite-size effects, and, as the chain length N — > 00, this maximum occurs at the CAP. 

Consider now the steepest descent treatment for the deep adsorption regime. With Eq. p2)) (after the rcscaling 
w S^w) in Eq. (PT|) . the saddle point becomes 
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where 62 = (1 - X)~'^^~^^ (recall that 1 - A = 7^=2 = 1.343). The main contribution comes from the exponential term 
in Eq. ([221) which is given by 



exp[-7Vlnz*«o) ^ {K+l)ln^o] = (^,^)A- ^^'^^''"'^li"-/)^!! (36) 



For A'^ — ii' ^ 1 the expression for the non-normalized PDF takes on the form 

, N-K 



^ iN-K)(^-^m-^) ^^^^ 



where bs is a constant of the order of unity. Normalization of this distribution yields 



" 7^(x)(iV-if)(l-^)(^-^) ^^^^ 
where the parameter x = b^ifJ-^l l^2w)N^^^ and the normalization constant 

^(X)- E (jV-j^)(i-A)(^-i.) - (39) 

The validity condition, Eg. ([23]) . in the deep adsorption regime (after substitution of Eq. (|32l) into Eq. (|23l) ) requires 

{N ~ Kf/^-^ > (^3//i2w)Ari"^, (40) 

i.e., K should not be very close to N . In Fig. [3^ we show the PDF of the number of chain contacts, P{K), for a 
free chain without pulling and several adsorption strengths of the substrate. One can readily verify that visually the 
shape of P{K) resembles very much a Gaussian distribution at moderate values of Ec ~ 1-7 < e < 6.0. The PDF 
variance goes through a sharp maximum at e > ec and then declines, as expected from Eq. p4)) . 

One should note that in the /-ensemble (where the force / and not the distance h acts as a controll parameter 
[1^ ) the order parameter n ~ K/N undergoes a jump at the detachment adsorption energy e_D. This means that 

7V(9^ lnz*/9e^)ej3 — —{dK/de)^^ — > —00. Thus at the detachment point the variance of the fluctuations {K — KY — 
—N{d'^ lnz*/9e^)co ~^ 00, which practically means that for chains of a finite length the distribution at e = en becomes 
very broad, in sharp contrast to Eq. ([M)) . This has indeed been observed in our MC-simulation results (see Fig. 12 in 
ref. M)- 



C. P{K) distribution in the subcritical regime w < Wc of underadsorption 

In the subcritical regime, w < Wc, the fraction of adsorbed points (order parameter) n = K/N = 0, in the 
thermodynamic limit. Nevertheless, K ^ and one can examine the form of the PDF P{K). Ai w <Wc the solution 
for z*{w) (the simple pole of S(z, w) in the complex z-plane) does not exist because V{wz)U{z) < 1 (see Ea. ([T5)) ). 
However, the tail GCE-partition function Q{z) — 1-|-<I>(1 — 71, ii^z) oc T (-ji) / {1 — fj.^z)'^^ has a branch point ed, z — l//i3 
(see Eq. (A 11) in ref. [I3|) which governs the coefficient at z^ , i.e., the partition function En{w). The calculation 
(following Section 2.4.3 in ref. [2^) yields 

1 - C(a) $(A,^2w/M3) 

where A = 1 — 7^=2 and we have also used that at z — l/fJ-3 the loop and train GCE-partition functions are 
[7(1/^3) = <i>(a, 1) — ({a) and ¥{11]/ fi^) = ^{X, ^2w/ ^j.^) respectively. This expression has a pole at w = Wc (cf. Eq. 
(^5)1 ) which yields the coefficient of , i.e. '^n,k- Recall that SAr(w) = Y^^=o'^n,kw^ , so that P{K) oc 'E.^.kw^ ■ 
Expansion of the denominator in Eq. (|4ip around w = Wc reveals the simple pole as follows 



-n{w) = 1— . (42) 

C(a) <I'(-7<i=2, M2Wc/m3) ywc--- 
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In (|42)) we have used the relationship x{d/dx)^{l — ^d=2,x) = ^{~jd=2,x). The coefficient of , i.e., ^n,k, 
is proportional to . Therefore P{K) cx {w/wc)^ — cxp[— (e^ — ^)K]. Taking the normalization condition 
^^^Q P{K) — 1 into account, the final expression for P{K) can be recast in the form 

P{K) = exp[-if(., - e)l (43) 

1 - exp[-iV (ec - e)J 

i.e., P{K) has a simple exponential form. The calculation of the average K — '}2!k=qKP{K) leads to the simple 
result K = [exp(ec — e) — 1]^^ ~ ^li^c — e), i-e., if — > oo at e ^ Cc- On the other hand we know that K = N'^ at 
e — Ec- In order to prevent a divergency at e = ec, one should incorporate an appropriate cutoff in the PDF given by 
Ea. p5)) . With this the distribution is given by 

-<«^T^^i^^|-^r^-l-«.-..i/»*)l^ (4.) 

Thus the expression for the average number of adsorbed monomer has the correct limit behavior, i.e., 



D. Probability distribution function P{K) in the ft-ensemble 



Eventually we examine how the fixed chain-end hight h affects the PDF of the number of contacts K. To this 
end we refer again to Fig. [2k where an adsorbed chain with a fixed height h of the last monomer is depicted. The 
adsorbed chain consists of a tail of length M and of an adsorbed part with N — M beads. One should bear in mind 
that M is function of the control parameters h and oi w = exp(e) given by eas.l|12p and (|13p. The partition function 
of the adsorbed part is then given by 

N 

■~n{w)^Y^ Sa,,kw^ =exp{-[iV-M(/i,w)]lnz*(w)} (46) 

A'=0 

where we took into account that the free energy of the adsorbed Dortion is eiven as F=h» = knTlN — M(h. w)] lnz*(u') 
( 




FIG. 3: (a) Probability Distribution Function of the number of chain contacts with the grafting plane, K, for different degrees 
of adhesion e and no pulling force / = as follows from eqs. (|28p and (|38p . The variance of P{K) is shown in the inset for 
different e. (b) The same as in (a) for e = 1.75 and different heights of the chain end h. The change of the variance {K'^) — (K)^ 
with h is shown in the inset. 



As mentioned above, the PDF P{K) oc :=.n,k , so that by means of rescaling w £,w and Wc — > £,Wc the PDF 
can be found as the coefficient of i.e., 

. J_ I exp{-[iV-M(. eH]ln.-(e^)} (4,) 
27ri / ^«+r 
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As before, the steepest descent method can be used to calculate the integral in Eq. (|47)) . However, in this case the 
calculations are more complicated and we have relegated most of them to Appendix B. As may be seen there, the 
saddle point equation can not be solved analytically in the general case but could be treated iteratively for relatively 
small heights h. 

One can readily see that for ft, = Eq. (|B7p reduces to Ea.lfSS]). The PDF, following from Eq. jBTj is shown in 
Fig. [3)3 for e — 1.75 and several values of the height h. It can be seen that the curve for ft = coincides with the 
curve for e — 1.75 in Fig. [3^ as it should be. Evidently, both the mean value K and K"^ — K decline with growing ft. 



IV. MONTE CARLO SIMULATION MODEL 

We use a coarse grained off-lattice bead-spring model [l^l which has proved rather efficient in a number of polymers 
studies so far. The system consists of a single polymer chain tethered at one end to a flat impenetrable structureless 
surface. The surface interaction is described by a square well potential, 

(48) 

The strength e is varied from 2.0 to 5.0 while the interaction range Tc = 0.125. The effective bonded interaction is 
described by the FENE (finitely extensible nonlinear elastic) potential: 




Ufene = -K{1 - lofln 



l-lo 

^max ^0 



(49) 



with K = 20, and the fully stretched-, mean-, and minimum bond lengths l„iax = 17^0 = 0.7, ^„ij„ = 0.4. The 
nonbonded interactions between monomers are described by the Morse potential: 

^'^'^^^^ = exp(-2Q:(r - r™„)) - 2exp(-Q;(r - r-rmn)) (50) 
Cm 

with a ~ 24, r„iin — 0.8, e]\j/kBT = 1. In few cases, needed to clarify the nature of the polymer chain resistance 
to stretching, we have taken the nonbonded interactions between monomers as purely repulsive by shifting the Morse 
potential upward by cm and removing its attractive branch for r > rmin- 

We employ periodic boundary conditions in the x — y directions and impenetrable walls in the z direction. The 
lengths of the studied polymer chains are typically 64, and 128. The size of the simulation box was chosen appropriately 
to the chain length, so for example, for a chain length of 128, the box size was 256 x 256 x 256 . All simulations 
were carried out for constant position of the last monomer z-coordinate, that is, in the fixed height ensemble. The 
the fluctuating force /, exerted on the last bead by the rest of the chain was measured and average over about 
2000 measurements. The standard Metropolis algorithm was employed to govern the moves with self avoidance 
automatically incorporated in the potentials. In each Monte Carlo update, a monomer was chosen at random and 
a random displacement attempted with Az, Ay, Az chosen uniformly from the interval —0.5 < Ax, Ay, Az < 0.5. 
If the last monomer was displaced in z direction, there was an energy cost of — /Az due to the pulling force. The 
transition probability for the attempted move was calculated from the change AC/ of the potential energies before 
and after the move was performed asW = exp{—AU/kBT). As in a standard Metropolis algorithm, the attempted 
move was accepted, if W exceeds a random number uniformly distributed in the interval [0, 1). 

As a rule, the polymer chains have been originally equilibrated in the MC method for a period of about 5 x 10^ 
MCS after which typically 500 measurement runs were performed, each of length 2 x 10^ MCS. The equilibration 
period and the length of the run were chosen according to the chain length and the values provided here are for the 
longest chain length. 



V. MONTE CARLO SIMULATION RESULTS 



In order to verify the theoretical predictions, outlined in Section|IIl we carried out extensive Monte Carlo simulations 
with the off-lattice model, defined in Section IIVI In these simulations we fix the end monomer of the polymer chain 
at height ft above the adsorbing surface, and measure the (fluctuating) force, needed to keep the last bead at distance 
ft, as well as the corresponding fraction of adsorbed monomers n. These computer experiments are performed at 
different strengths e of the adsorption potential, Eq. (HS)) . In Figs. 2^, b we compare the predicted dependence of the 
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FIG. 4: (a) Order parameter (fraction of adsorbed monomers) n variation with changing height h/loN of the fixed chain-end 
for polymers of length = 128 and different adsorption strength e/kBT. (b) Variation of n with e/kBT for different fixed 
positions of the chain-end h/loN as it is seen from MC-data. Insets show the resulting n — e relationship at several fixed heights 
h. 



order parameter n on the (dimensionless) height h/l^N at several values of 2.0 < e < 5.0 with the results from MC 
simulations. Note, that the critical point of adsorption Cc ~ 1-69 so we take our measurements above the region of 
critical adsorption. Typically, both in the analytic results. Fig. 2^, and in the MC-data, Fig. Uja, for N = 128, one 
recovers the predicted linear decrease of n with growing h. Finite-size effects lead to some rounding of the simulation 
data (in Fig. [4|3 these effects are seen to be larger for TV = 64 than for N = 128) when n ^ so that the height of 
detachment /id is determined from the intersection of the tangent to n{h) and the a;— axis where n = 0. Evidently, 




FIG. 5; Phase diagram showing the dependence of the critical height of polymer detachment from the substrate, ho/loN, with 
the relative strength of adsorption [e — Sc) /ksT , where SclksT is the critical point of adsorption at zero force. The theoretical 
curves follow Eq. (fT7|l . 

with growing adsorption strength, e, larger height ho is needed to detach the polymer from the substrate. Thus, one 
may construct a phase diagram for the desorption transition, which we show in Fig. [S] The theoretical prediction is 
given by eq. (fT7|) . 

In the insets of Fig.[3^,b we also show the variation of the fraction of adsorbed segments n with adsorption strengths 
e for several heights 20 < /i < 50 of a chain with N — 128. It is evident that, apart from the rounding of the MC data 
at n — > 0, one finds again good agreement between the behavior, predicted by Eq. (fTS)) . and the simulation results. 

The gradual change of n in the whole interval of possible variation of h suggests a pseudo-continuous phase transition, 
as pointed out in the end of Section [Til Eq. (fTB)) . Of course, if h is itself expressed in terms of the measured pulling 
force, one would again find that n changes abruptly with varying / at some threshold value /d, indicating a first 
order transformation from an adsorbed into a desorbed state of the polymer chain. 
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— e=1.3 




FIG. 6: PDF of the order parameter (fraction of contacts with the plane) for different adsorption strength e at zero force. The 
critical adsorption point (CAP) Ec ~ 1.67. The change of the variance (n^) — (n)^ with varying e is displayed in the inset. 

It has been pointed out earlier by Skvortsov et al.[l3] that, while both the fixed- force and the fixed-height ensembles 
are equivalent as far as the mean values of observables such as the fraction of adsorbed monomers and other related 
quantities are concerned, this does not apply to some more detailed properties like those involving fluctuations. 
Therefore, it is interesting to examine the fluctuations of the order parameter, n, for different values of our control 
parameter h, and compare them to theoretical predictions for P{K) from Section [ill First we compare the order 
parameter distribution P{n) for zero force. Fig. [HI obtained from our computer experiment, to that, predicted by 
Eqs. ((28| . ([38)) . (|44p. and displayed in Fig. [Sji. It is evident from Fig. [6] that for free chains at different strengths 
of adhesion there is a perfect agreement between analytical and simulational results. For rather weak adsorption 
e = 1.3 ± 1.6 < Ec = 1.67 in the subcritical regime, one can verify from Fig. [6] that P{n) gradually transforms from 
nearly Gaussian into exponential distribution, as expected from Eq. (|^ . For e > Ec the distribution width ((n — (n))^) 
grows and goes through a sharp maximum in the vicinity of Cc, and then drops as e increases further - compare insets 
in Fig. [6] and Fig. EK. 




FIG. 7: (a) Probability distribution P{n) of the order parameter n (i.e., the fraction of adsorbed monomers) for A'' = 128 
and e = 3.0 at different heights of the chain-end h over the grafting plane. In the inset we show P{n) at the detachment line 
ho ~ 54.3. (b) Variation of the second- and third central moments of P{n) with h. The maximum of ((n — {n))'^) is reached 
at /i = ho- 

Let us consider now PDF in the presence of pulling. In Fig. [7^ we display the distribution P{n) measured in the 
MC simulations for different heights h and constant adsorption energy e = 3.0. One can readily verify from our 
results that far enough from the detachment line, h < hn, the shape of P{n) looks like Gaussian and that the second 
moment, ((An)^), remains unchanged with varying height h. Of course, when h hjj, the maximum of P{n) shifts 
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to lower values of n. Only in the immediate vicinity of ha, where n ^ and the fluctuations strongly decrease, 
one observes a significant deviation from the Gaussian shape - cf. the inset in Fig. [7^. The latter is illustrated in 
more detail in Fig. [TId where we show the measured variation of the second moment, ((An)^), and that of the third 
moment, ((n — (n))^) with increasing height h. The deviation from Guassinity in P{n), measured by the deviation of 
the third moment from zero, is localized in the vicinity of the detachment height ho- The corresponding theoretical 
prediction for the relatively weak adsorption strength is depicted in Fig. [Sb. It can be seen that with increasing h 
the almost Gaussian distribution tends to Poisson-like one. Also the fluctuations decrease with h in accordance with 
MC-findings. 

The force /, exerted by the chain on the end-monomer, when the latter is kept at height h above the surface, is 
one of the main properties which can be measured in experiments carried out within the fixed-height ensemble. Note 
that / has the same magnitude and opposite sign, regarding the force, applied by the experimentalist. The variation 
of the force / with increasing height h is shown in Fig.[5|i for several values of the adsorption potential 2.0 < e < 5.0. 
In Fig. [8^ we distinguish between two contributions to the total force /, acting on the end bead. The first stems from 
the quasi-elastic forces of the bonded interaction (FENE) whereas the second contribution is due to the short-range 
(attractive) interactions between non-bonded monomers (in our model - the Morse potential) . A typical feature of the 
f — h relationship, namely, the existence of a broad interval of heights h where the force remains constant (a plateau 
in the force) is readily seen in Fig. [8^. With growing strength of adsorption e the length of this plateau as well as the 
magnitude of the plateau force increase. Note, that for e = no plateau whatsoever is found. Upon further extension 
(by increasing h) of the chain, the plateau ends and the measured force starts to grow rapidly in magnitude - an 
effect, caused by a change of the chain conformation itself in the entirely desorbed state. 

A closer inspection of Fig. [5^ reveals that the non-bonded contribution to /, which is generally much weaker than the 
bonded one, behaves differently, depending on whether the forces between non-nearest neighbors along the backbone 
of the chain are purely repulsive, or contain an attractive branch. While for strong adsorption, e > 3.0, a plateau 
is observed even for attractive non-bonded interactions, for weak adsorption, e < 2.0, an increase of the non-bonded 
contribution at h/l^N w 0.35, (seen as a minimum in Fig.jSK) is observed. This effect is entirely missing in the case of 
purely repulsive nonbonded interactions - see the inset in Fig. [5^ where the contributions from bonded and non-bonded 




FIG. 8: (a) Variation of the two components to the total force, exerted by the chain on the end-monomer which is fixed 
at (dimensionless) height h/loN for different adsorption potentials 2.0 < e/ksT < 5.0: bonding interactions (full symbols) 
and non-bonding Morse interactions (empty symbols). In the inset the same is shown for a neutral plane e — 0.0 for purely 
repulsive monomers (triangles) and for such with a weak Morse attraction (circles), (b) Variation of the total force (plateau 
hight) exerted by the AFM tip on the chain-end for chain length N = 128 with the relative strength of adsorption [e ~ Ec) /ksT . 

interactions are shown for a neutral surface e = 0. If one plots the magnitude of the measured force at the plateau 
against the corresponding value of the the adsorption potential, e, one may check the theoretical result, Eq. ([T2|) - 
Fig.[Hj3. Evidently, the theoretical predictions about /p agree well with the observed variation of the detachment force 
in the MC simulation, both within the BS- or FJBV models, as long as only the excluded volume interactions in the 
MC data are taken into account. If the total contribution to /p, including also attractive non-bonded interactions in 
the chain, is depicted - black triangles in Fig. [Sh - then the agreement with the theoretical curves deteriorates since 
the latter do not take into account the possible presence of attractive non-bonded interactions. 

The f — h relationship, which gives the equation of state of the stretched polymer, may be derived within one of the 
different theoretical models, e.g., that of BS-, Eq. ([H]), or FJBV-model, Eq. ([5]), as mentioned in Section|lTl Which of 
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FIG. 9: (a) Variation of the total applied force / with growing height of the end monomer in terms of Pincus reduced variables, 
floN" /ksT versus h/loN" , for a polymer with purely repulsive nonbonded forces for A'^ = 64, 128. (b) The same as in (a) but 
in terms of reduced units flo/ksT versus h/l^N for purely repulsive (empty symbols) as well as for usual Morse potential (full 
symbols) of nonbonded interactions between monomers. The FJBV- model results, Eq. ([Sj, is shown by a solid line. Arrows 
indicate the unperturbed gyration radius positions Rg/N for A'' = 64, 128 

these theoretical descriptions is the more adequate can be decided by comparison with experiment. In Fig. [9^,b, we 
present such comparison by plotting our simulation data using different normalization for the height h. From Fig. [5^ 
it becomes evident that the data from our computer experiment for iV = 64 and iV = 128 collapse on a single curve, 
albeit this collapse only holds as long as h/loN'^ < 3.0 for the BS-model while it fails for stronger stretching. In 
contrast, this collapse works well for all values of h, provided the height is scaled with the contour length of the chain 
A'', rather than with N'^, as in the FJBV model - Fig.[S}3 - regardless of whether a purely repulsive, or the full Morse 
potential (which includes also an attractive part) of interactions is involved. The analytical expression, Eq. ([5]), is 
found to provide perfect agreement with the simulation data for strong stretching, h/l^N > 0.4. From the simulation 
data on Fig. [5] one may even verify that the force / goes through zero at some height h > and then turns negative, 
provided one keeps the chain end very close to the grafting surface (cf. eq.(l6|)). 

VI. SUMMARY 

In the present work we have treated the force-induced desorption of a self-avoiding polymer chain from a flat struc- 
tureless substrate both theoretically and by means of Monte Carlo simulation within the constant-height ensemble. 
The motivation for this investigation has been the necessity to distinguish between results obtained in this ensemble 
and results, derived in the constant-force ensemble, considered recently fl^, as far as both ensembles could in principle 
be used by experimentalists. We demonstrate that the observed behavior of the main quantity of interest, namely, the 
fraction of adsorbed beads n (i.e., the order parameter of the phase transition) with changing height h differs quali- 
tatively from the variation of the order parameter when the pulling force is varied. In the constant-height ensemble 
one observes a steady variation of n with changing h whereas in the constant-force ensemble one sees an abrupt jump 
of n at a particular value of fu, termed a detachment force. However, this should not cast doubts on the genuine 
first-order nature of the phase transition which can be recovered within the constant-height ensemble too, provided 
one expresses the control parameter h in terms of the average force /. This equivalence has been studied extensively 
for Gaussian chains by Skvortsov et al. fT3| who noted that ensemble equivalence does not apply to fluctuations of 
the pertinent quantities too. 

Indeed, in our earlier study we found diverging variance of the PDF P{n) at fjj whereas in our present study 
the fluctuations of the order parameter are observed to stay finite at the transition height hu. These findings confirm 
theoretical predictions based on analytic results which we derive within the GCE-approach. Within this approach we 
have explored two different theoretical models for the basic force - extension relationship, namely, the bead-spring 
(BS) model as well as that of a Freely- Jointed Bond- Vectors (FJBV) model. Our simulation results indicate a good 
agreement between theory and computer experiment. 
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APPENDIX A: FREELY JOINTED BOND VECTORS MODEL 



The deformation law in the overstretched regime (when the chain deformation is close to its saturation) could be 
treated better within the FJBV model. Consider a tethered chain of length N with one end anchored at the origin of 
the coordinates and an external force /at acting on the free end of the chain. The corresponding deformation energy 
reads 



N 



-Jn rj^ = -fN ^ 6iCos( 



(Al) 



where rj^ is the z-coordinate (directed perpendicular to the surf ace) of the chain end, bi and 9i are the length and 
the polar angle of the i-th bond vector respectively. The corresponding partition function of the FJBV model is given 
by 
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The average orientation of the i-th bond vector can be calculated as 



< cos Hi >- 
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From Eq. (|A3|) the chain end mean distance from the surface, h, is given by 



h^J2b^< cose, bNC (^^^ 



(A2) 



(A3) 



(A4) 



where we have taken into account that the lengths of all bond vectors are of equal length, hi — h, and C{x) 
coth(a;) — 1/x is the Langevin function. This leads to the force - distance relationship 



N 



kRT 



h 

m 



, at h/bN < 1 
, at h/bN < 1 



2h 
bN 



(A5) 



which we use in Sec. II. The notation C ^{x) stands for the inverse Langevin hmction. 



APPENDIX B: CALCULATION OF PDF IN THE ft-ENSEMBLE 



Using Eg. ([24]) for In z*{w) as well as Eqs. (IT2]) and (for the BS-model), one obtains an expression for the tail 
length 



af{w - Wc) ^ 

The saddle point (SP) equation in this case reads (cf. Eg. pT]) ) 



K + 1 ... . y d\nz* 



dM 



(Bl) 



(B2) 
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Taking into account Eqs. (|24|) and (jBip . after introducing the notation y ~ {w — Wc)^q, the SP-equation can be recast 
into 



h 



h \ In /13 



IqN J y(i-'')/</' 



(B3) 



where Bi and i?2 are constants of the order of unity. In the particular case h = Eq.jBSj) goes back, as expected, to 
Eq.lUe]). Eq. ((B3)) can be solved iteratively for h < IqN as 
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KY f h 
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(B4) 



As before, the main contribution in the integral given by Eq. (|47p reads 
P{K) oc cxp{-[N - M{h,^ow)]\nz*{^ow) " {K + l)ln^o} 



{w — Wc)N'^ exp 



1 - 



(h/loN) 



{K/Ny-[piK/N, h/loN)] 



K 



piK/N, h/loN) \ [Kp{K/N, h/loN)]-'^\B5) 



where one introduces the notation 



p{K/N,h/l,N) = l + Bij^jj^ + B2\np3 ' 



{K/Ny-'^ 



(B6) 



After normalization, the final expression for the PDF reads 



P{K) 



W{v,h/loN) 



1] exp 



{h/loN) 
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p{K/N, h/loN) } [Kp(K/N, h/loN)]-^''{B7) 



where rj = {w — Wc)N'^ as before and the normalization constant reads 
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W{rj, h/loN) =Y, S ^ exp 
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p{K/N, h/loN) } [Kp{K/N, h/loN)]-^'' (B8) 



Again one can readily see that for h = Eq. (|B7|) reduces to Eq.([28|). The PDF, following from Eq. (jB7| is shown in 
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Fig. [3)3 for several values of the height h. Evidently, both the mean value K and the variance K-^ — K decline with 
growing h. 



[1] S. B. Smith, Y. Cui, and C. Bustamante, Science 271, 795(1996). 

[2] T. J. Senden, J.-M. di Meglio, and P. Auroy, Europ. Phys. J. B3, 211(1998). 

[3] T. Hugel, M. Grosholz, H. Clausen- Schaumann, A. Pfau, H. E. Gaub, and M. Seitz, Macromolecules, 34, 1039(2001) 
[4] H. Hugel, M. Seitz, Macromol. Rapid. Commun. 22, 989 (2001). 
[5] A. Rohrbach, E.H.K. Stelzer, J. Appl. Phys. 91, 5474 (2002). 

[6] L. Sonnenberg, Y. Luo, H. Schaad, M. Seitz, H. Golfen, H.E. Gaub, J. Am. Chem. Soc. 129, 15364 (2007). 
[7] C. Friedsam, A. Del Campo Becares, M. Seitz, and H. E. Gaub, New J. Phys. 6, 9(2004). 
[8] S.K. Kufer, E.M. Puchner, H. Gumpp, T. Liedl, H.E. Gaub, Science, 319, 594 (2008). 
[9] M. Rief, F. Oesterheh, B. Heymann, and H. E. Gaub, Science 275, 1295(1997). 
[10] C. Ortiz and G. Hadziioannou, Macromolecules 32, 780(1990). 

[11] M. Grandbois, M. Beyer, M. Rief, H. Clausen-Schaumann, and H. E. Gaub, Science 283, 1727(1999). 

[12] A.M. Skvortsov, L.I. Klushin, T.M. Birshtein, Polymer Sci. (Moscow) Ser. A, 51, 1 (2009). 

[13] S. Bhattacharya, V.G. Rostiashvih, A. Milchev, T.A. Vilgis, Phys. Rev. E 79, 030802(R) (2009). 

[14] S. Bhattacharya, V.G. Rostiashvih, A. Milchev, T.A. Vilgis, Macromolecules, 42, 2236 (2009). 

[15] C. Vanderzande, Lattice Model of Polymers, Cambridge University Press, Cambridge, 1998. 

[16] E. Eisenriegler, K. Kremer, K. Binder, J. Chem. Phys. 77, 6296 (1982). 

[17] E. Eisenriegler, Polymers Near Surfaces, World Scientific, 1993. 

[18] T. Kreer, S. Metzger, M. Miiller, K. Binder, J. Chem. Phys. 120, 4012 (2004). 

[19] A.Yu. Grosberg, A.R. Khokhlov, Statistical Physics of Macromolecules (AIP, New York, 1994). 



17 



[20] J. dos Cloiscaus, G. Jannink, Polymers in Solution, Clarendon Press, Oxford, 1990. 
[21] P.-G. dc Germes, Scaling Consept in Polymer Physics, Cornell Univercity Press, 1979. 
[22] Y.-,]. Sheng, P.-Y. Lai, Pliys. Rev. E 56. 1900 (1997). 

[23] M. Wittkop, J.-U. Sommer, S. Kreitmeier, D. Goritz, Phys. Rev. E 49, 5472 (1994). 
[24] J.U. Schurr, S.B. Smith, Biopolymers, 29, 1161 (1990). 

[25] J. Rudnick, G. Gaspari, Elements of Random Walk, Cambridge Univercity Press, Cambridge, 2004. 

[26] S. Bhattacharya, H.-P. Hsu, V.G. Rostiashvili, A. Milchev, T.A. Vilgis, Macromolecules, 41, 2920 (2008). 

[27] K. Binder and A. Milchev, J. Computer-Aided Material Design, 9, 33(2002) 



